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~qq ". We investigate nonperturbative dilatonic solutions of the wide class of the modified 

gravity models including the Gauss-Bonnet terms with a general F(G) Lagrangian. 

CN ■ We show that presence of the Liouville-like solutions is a characteristic feature of 

these models. 
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m ■ 
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The discovery that the Universe is presently undergoing accelerated expansion [l| has in- 
spired several researchers to search for physical models which could explain this phenomenon. 
Among different theories applied for this purpose, one of the most interesting suggestions, 
beside of use of the cosmological constant [2] and quintessence [3J, is the idea of a modified 
gravity according to which the common Hilbert-Einstein gravity action is extended by addi- 
Q\ ; tive terms depending on some scalars constructed on the base of the metric tensor. In the 
mostly used form, such a modification implies in the theory called F(R) gravity where the 
Lagrangian is a function of the scalar curvature (see e.g. [4j). The cosmological implication 
of such models were discussed, for instance, in j^J, and in our recent works ja, 0] some static 
OO ! solutions for such models were found. 

However, the F(R) modification is not the unique possibility to change the action of gravity 
without introducing extra fields. One of the popular suggestions for the modified gravity 
action, beside of the F(R) models and of the Weyl gravity is the Gauss-Bonnet gravity 
whose different classical aspects (including properties of some solutions of the equations of 
motion) were studied in a series of papers [9|. Recently, the Gauss-Bonnet gravity models 



were applied also for the brane studies [lOj, including holography aspects 



111 ], and solving 



the problem of cosmic acceleration, within which the modified Gauss-Bonnet gravity was 
treated as an alternative explication of dark energy . Also, the Gauss-Bonnet gravity 



models were applied to studies of higher- dimensional gravity |14j and of black holes 15(. Some 
observational estimations for Gauss-Bonnet gravity are given in [l7. \l8\. 

Another direction of current interest is related to cosmology. Since we will be working with 
a single space coordinate, our investigation is somehow similar to the study which generalizes 
the FRW scenario by including the Gauss-Bonnet contribution, since there one requires time 
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evolution, and time is similar to a _space coordinate, unless for its signature. Recent investi- 
gations on this issue appeared in 16(, and there one can also find the Liouville-like behavior 
which will appear in the present study. The string-inspired investigations done in 16j and 
in references therein further stimulates the present study, bridging its contents to cosmology 
and the up to now misterious dark energy fuel. 

The interest to modified gravity models naturally stimulates the search the interest for the 
exact solutions in these models. In our papers [s, 0] we have shown that the Liouville-like 
solutions, being the simplest examples of the exact solutions and, in certain cases, correspond- 
ing to the anti-de-Sitter space, are the typical one for the F(R) gravity models. Therefore, 
the very natural question is whether such solutions are possible in the more general modified 
gravity models, that is, in different versions of the Gauss-Bonnet-like models for the gravity 
in arbitrary space-time dimensions. We note that this interest is further supported by the fact 
that the time-dependent Liouville-like solutions were shown to arise also in the cosmological 
context fl6| . 

First of all, let us again, just as in [fjj], restrict ourselves to the simplest conformally flat 



space 



19] , with the dynamics is now concentrated in the conformal sector of the gravity, and 



the metric tensor g^x) is given by 



6 ^ "Htmi 



(1) 



with <j(x) standing for the conformal factor, the dilaton field. We use metric (+, — , . 
and work with dimensionless fields and coordinates, for simplicity. 
In this case, one can find that the Christoffel symbols look like 



pu 



(2) 



The Riemann curvature tensor is equal to 



R 



p.paf3 



pi "p f) V -I- F Y u F Y v 

u a L fi,pf3 u f3 L p,pa T 1 p,ua L pR 1 p,,uf3 L p, 



pa 



= e a [-(d a d p ar]^ - dpdpar]^) + -(dfid^arjap - d a dp,(jr]p p ) + 

1 1 

+ -{dpdpar] a p - d p d a cn]pp) + -{dpd a ar]p p - dpdporiup) + 

+ ^d x ad x a(r] a pr] l3 p - VfipVap)}, 

and this implies that the Ricci tensor in D-dimensional space-time is 

D-2 



Rpa Rpap 



[D - 2)d p d a a + r!p a Da 
Thus, the scalar curvature has the form 



d p od a o - rj pa d^adp(T 



(3) 



(4) 



R = 9 aP Rap = {D- l)e" CT 



Da ( ^-j^ ) iPaO^a 



(5) 



The useful contractions are also 

R.uXpR^" = e~ 2 °(D-2) 

+ d a ad a aDa - d a ad b ad a d b a] , 



d a d b ad a d b a + (Da) 2 + i(£> - l)(d a ad a a) 2 + 

o 



(6) 



and 



-2a 



+ h D -l){D-2)\cfad a af+ 
lo 

The generic invariant action, which looks like 

S = J d D x^\(AR 2a + B{R^ p R^) b + C{R^) C ) , 
in the dilatonic sector takes the following form 



(7) 



(8) 



1 



+ Be- 2ba 



, d a ad a a 
4 1 



2a 



+ 



(D - 2)d a d b ad a d b a + (Da) 2 + -(D-1)(D- 2)(d a ad a a) 2 + 

o 

+ (D - 2)d a ad a aDa - (D - 2)d a ad b ad a d b a] 26 + 

+ Ce- 2ca 



( D — 2) 2 3D 1 

1 4 ) d a d b ad a d b a + (— - 1)(D^) 2 + Yq( D ~ 1)P - 2) 2 (d a ad a a) 2 + 



+ (D - 2)(2D - Z) daadaaDa _ (D^T^rto 



(9) 



The important case of this generic action is the Gauss-Bonnet action depending on the Gauss- 
Bonnet invariant of the form 



G — R 2 — AR^ V R^ U + R^ v \ p R^ uXp ', 
with the explicit expression for it having the form 



(10) 



G = e- 2a {D -2)(D -3) 



[na) 2 - d a d b ad a d b a + 



(jD _1 )(jD _ 4) 
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(d a ad a a) 2 + 



+ d a ad b ad a d b a + ^ ^d a ad a aDa 



(11) 



In particular, in four dimensions, the Gauss-Bonnet action is a surface term: 

S GB = J d 4 x^\G = J d 4 xe 2a G = J d 4 xd a (2(dVQr - d a d b ad b a) + d a ad b ad b a) . (12) 
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In two and three dimensions, the Gauss-Bonnet invariant is an identical zero. 

In this paper, we will concentrate on the domain wall case. To do this, we suggest that the 
dilaton a depends only on one spacial coordinate, say z. So, (Oa) 2 = (u") 2 , d a db<rd a d b cr = 
(a") 2 , {d a ad a af = (a') 4 , d a ad b ad a d b a = (a') 2 a", 8 a ad a aUa = (a')V'. 

First of all, we can show that the Liouville-like kink solutions are characteristic ones for a 
wide class of dilaton gravity models, with we suggest the restriction up to the domain wall 
dynamics. Indeed, in this case the scalar curvature is given by, 



R 



D 



a 



-a 



and the two other scalars look like 

R»vX P R^ Xp = (D - \)e~ 2a 



o") 2 + UD-2){e') A 



and 



R^BT — (D — l)e 



-2(7 



D 

7 
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(13) 



(14) 



(15) 



Let us consider the generic action (jHJ). This action, in the dilatonic sector, by use of the 
expressions above, can be expressed as 



S 



d D xe^ 

2ba 



Ae~ 2a ° ( ai (a") 2 + a 2 {a'Y + a 3 a"(a') 2 ) a + 

Be~™ (f3 1 (a") 2 + (3 2 (a') A + f3 3 a"(a') 2 ) b + 
+ Ce~ 2c ° (71 (a") 2 + 72(<x? + ^"WfT 



(16) 



where ji, A, B, C are some constants whose exact form is not important. 

We can obtain the equations of motion for this theory. We note that inside each paren- 
theses, the number of derivatives in each term is equal to 4, thus, the number of derivatives 
in any term in the equations of motion proportional to A is 4a, to B is 46, to C is 4c. Thus, 
we find the equations of motion in the form 

Ae (-2a + f } , a klk2k3 n (arvr 2 (o fc3 + 

ki,k2,k3 fci+2fc2+3fc3=4a 



+ 



+ Ce(~ 2c+ TV c mi 



"'2 m 3 W ' (1 

mi,m2,m-j mi+2m2+3m3=4c 



') rai ((iT 2 (ff'T 3 = o, 



;i7) 



where a kl k 2 k 3 ^ b niri2n3 , c mim2m3 are some numerical coefficients. Then, we introduce the known 
Liouville-like ansatz 



a' = Ne la , 



(18) 
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which implies in 

a" = N 2 le 2lcT , a'" = 2N 3 l 2 e 3l<7 . (19) 

Substituting these ansatzes, we see that in the term proportional to A, the field dependent 
multiplier factorizes giving the factor e (( 4/ - 2 ) a +T) cr ) whereas in the terms proportional to B 
and C we get, respectively, e^ Al ^ b+ ~^ a and e^ 4 '~ 2 - )c+ '?- )cr . As a result, the equation of motion 
is reduced to the algebraic one 

Ae (Al-2)a, N . a J- a klk2k3 J] l k ^ + 

ki,k2,ks fei+2fc2+3fc3=4a 

ni,ri2,n3 ni+2n2+3ri3=4b 

+ Ce^ l - 2 ^N 4c c ™™ J] r2+2m3 = °- ( 2 °) 

mi,m2,m3 mi+2m2+3m3=4c 

The different solutions of this algebraic equation for the coefficients can be found. The 
consistent cases are a = b = c, or I = |. Also the consistent cases are when any two of the 
three coefficients A, B, C are equal to zero, or one of them is equal to zero whereas the two 
others carry the same exponential factors, e.g., A = and b = c. 

The case I — 1/2 is of special importance. Indeed, in four dimensions the scalar curvature 
(jl"3l) for this case looks like 

R = -3iV 2 , (21) 

which is a negative constant. Thus, in four dimensions the case 1 — 1/2 corresponds to anti-de 
Sitter space, with the coefficients of (120]) are now related. 

We can also focus on another model. Let us find and solve the equations of motion for the 



model with a bit different structure, that is, with the action initially introduced in 13j (see 
also [12j, [20] for study of some aspects of this theory) : 

S = J d A x^W\{R + f{G)). (22) 

n 

Here we restrict ourselves by four dimensions. For simplicity, let us first, just as in [6|, consider 
the domain wall case. Thus, we get G = e" 2f7 ((cr') 3 ) / and R = — |e _fT (2cr" + (cr') 2 ) ^ s a first 
attempt, we study the action 



S = J d 4 x^/\i\(aR + bG 2 ) ~ J d 4 x 



96e- 2 V') V) 4 - + (o'f) 



(23) 



Here the sign ~ denotes that we restrict the coordinate dependence to the domain wall case. 
The corresponding equations of motion look like 

6&e- 2< V') V) 4 " Qb{{a'fa"e- 2a )" + 12b((a") 2 (a') 3 e' 2r7 )' + 
+ ae^a" + lae^a') 2 = 0. (24) 
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The most natural ansatz is (ITS]) implying in (|T9l and in 



(ct") 2 - a"(a') 2 = 1(1 - l)N 4 e 4la , (e mCT )' = Nme {m+l)c 
(e ma )" = N 2 m(m + l)e {m+2l)a . 



(25) 



Its substitution yields 

12 /WV 6 e (6/ - 3),7 [iV(7Z - 2) (6/ - 2) - l(Ul - 3)] - aN 2 a {21 + 1) = 0. (26) 

We have some distinct possibilities: first, we can choose 6 = 0, which yields I = —1/2, leading 
to usual gravitation; second, we can choose a = 0, which yields N = jjj^^zfj', and third, we 
can choose / = 1/2 with a, b ^ 0, so the exponential factor goes to 1, that is, we reproduce the 
adS solution, and we have the quintic algebraic equation, 3bN A (3N — 4) = 2a, which allows 
to find iV once a and b are given. 



Our next step is the model of the form (see f.e. 12]): 



S = J d 4 x^\{aR n + bG 



which generalizes the model (|22|) . In the domain wall case, it can be presented as 



S = J d 4 xe 2 ° ^[e- CT (-2a" - (ct') 2 )]™ + o[e~ 2 V) °"\ 
with a = (|) n a,b = 3 m 6. In this case, the equation of motion looks like 



{2 - n)e {2 - n) ° {-2a" - {a') 2 ) n + 



+ 2n[- e °( e -^>(-2o" - (ct') 2 )"- 1 )" - (e^ 1 ) \-2a" - [a') 2 )^ 1 )' eV] 



(27) 



(28) 



+ b(m - l)e {2 ~ 2m)u ( -{a') 2 a" 



1 \ m ~ ^ f 1 

■i- / #\ o #» \ I 1 / / i\ o n\ 3 , TTJ / / \ 3 / //\ 2 



4 ((-)V') d + y(a') 3 K) 2 



^ ^^4/ 7/\2 // //\2 | I i m ( m 2) / 2 

-(a ) (tr ) ((a ) + a a J + (a ) 



4 



<j I + oct a a H — ( ct J ct 



4 



( ) + 2^ ) CT 



;ct') 4 (ct") 2 [(ct") 2 + ct'ct"'] 



0. 



(29) 



The case a = we will consider in details further, the case 6 = can be read off from [6j, so 
we are interested only in the case when both a and b differ from zero. To solve this equation, 
we carry out the substitution used in fll8ll25p . Also, we introduce the notations 



T 
f 



16 2 ' 



D-A 
16 



-:v 4 + -iv 3 



(30) 
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Thus we get: 

~ e (2-„+2in)- Ar 2n ^ 2 _ n ^_ 2 i _ i) + 2n(-2Z - l)"" 1 ^ - l)(n - l)(2/n - / - n + 2 



(2 - 2m)T - 4/m(m - 1)T - 12mf+ 



m 



+ ^-N 2 l 2 {m - l)[16ra - 32 + 20Z 2 ] + (8m 2 - hm)N 2 l 2 



0. 



(31) 



It is easy to see, that if n = 2m, any / is possible. Otherwise, there exists an adS solution 
I = ~, for which any m and n are possible. Other parameters can be found by solving the 
algebraic equations. 

The natural continuation is a possible domain wall solution for the pure Gauss-Bonnet 
action in an arbitrary space-time dimension. In the domain wall case, the Gauss-Bonnet term 
is reduced to 

16 K J 2 K 

Let us now consider the case D > A (which is important from the extra dimensions viewpoint), 
and study the model with the action 



G = e- 2a (D -!)(£> -2)(D -3) 



(32) 



S 



d D xJ\i\G n = / d D xe^- 2n)a 



16 K J 2 K J 



(33) 



where we have omitted the irrelevant overall factor. The equations of motion look like 



^^(-') 4 + ^')V' 



n— 3 



D 



2n) 



— n[n — 1) 
D — A, 



x 



4 



16 V ; 2 y ' 



a ) a + a [a ) + -[a ) a 



(T') d + <7V 



X 



(34) 



16 v ; 2 V ; 



(a) a + 



n{n- l)(w-2) , a 
+ - o" 



D — 4 1 
—^—(a)(j +(J{a) +-{cr)a 



+ 



x 



nfn — 1"! 



A2 



16 V ; 2 V ; 



; cM a H — a a 



+ 2n(n - 1) 



(D-A). ,. 4 1 

16 (^) + 2^ r<r " 



x 



r ll\6 \ n I H Ht \ t l\l h 
[a J +3(7 a cr + — (cr ) cr 



4 (<T ) Cr + Cr ((T ) + -((7 J (7 



/ // 
cr cr 



0. 



Again, we carry out the substitution given by ( TT8lf25T) . After introducing the notations (J30 
we find that this substitution gives the following condition: either T = 0, or 

,D v _ ~„, ,~ , n 



2n)T - AnT[l(n - 1)T + 3] + -N 2 l 2 (n - 1) [16n - 32 + 20/ 2 ] + 
8n 2 - 5n)iV 2 Z 2 = 0. 



(35) 
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From here, we can relate C and I, thus, we again have the Liouville-like kink solutions. 
Also, we can study the model with the action 

S=- J d D x^/\f\(R-aG), (36) 

whose explicit form in the domain wall case is 

S = (D-1)J d D x[e^ 2 -^[a" + ^^(a') 2 } + 

+ a(D -2)(D- 3)e^/ 2 ^[^-V) 4 + \{a') V']l . (37) 

16 2 J 

For D > 4 the corresponding equation of motion are 

a" + ^yV) 2 - y Q {D - 3)(D - 4)[5(L> - 4)(a') 4 + 48(<x') V] = 0. (38) 
Implementing again the Liouville-like ansatz ({TBI we get the following form of this equation: 

N 2 e 2ltT (D + 21 - 2) - — (D - 3)(D - 4)e (4 '- 1),7 iV 4 [5( J D - 4) + 48/] = 0. (39) 

16 

The possible solutions are: first, a = (that lead us to usual Einstein-Hilbert gravity), and 
second, the case I = 1/2, whose solution is 



^^^S^) 1 (40) 



with 



N ^-^ ai D- 3 )(D-4) < 41 » 

Thus, we again have found the Liouville-like solutions. 

In summary, in this work we have studied a very wide class of the gravity models whose 
action depends not only on scalar curvature but also on the Gauss-Bonnet invariant and on 
other invariants constructed on the base of the Riemann curvature tensor, being considered 
both in four-dimensional space-time and in higher dimensions. As we have shown, all the 
models admit Liouville-like static kink solutions. Therefore, we can conclude that the presence 
of such solutions is common for different gravity models. Moreover, we succeeded to show 
that all these models in four-dimensional space-time allow the adS solution, and this seems 
to indicate that the systems admit supersymmetric extensions. 

We have considered solutions of the gravitational field equations which depend on a single 
spacial coordinate. It is natural to expect that in cosmology, where the metrics (and therefore 
the dilaton) depends only on time, would display similar behavior implying in the appearance 
of Liouville-like solutions. Within this context, we notice that the time-dependent Liouville- 



like solutions were shown to arise in the gravity context earlier, for instance, in 161 ]. a fact 
that extends the interest of this work to the cosmological scenario. 
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A very natural continuation of the present study could be the coupling of the dilaton 
to extra scalar matter fields. Such an extension could be applied to study more realistic 
scenarios, for example, in the braneworld and in the brane/cosmology contexts, giving thus 
a natural generalization of [7| . 
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